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Abstract 

We study the effects of next-to-leading order corrections on the evolution of the twist-two 
non-forward parton distribution functions in the flavour non-singlet sector. It is found that the 
deviation from leading order evolution is small for all values of the parton momentum fraction 
variable for moderately large values of the scale parameter. 
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1. Introduction. Recently there was significant progress in the perturbative QCD approach to 
deeply virtual Compton scattering — a process which allows to access the so-called non-forward 
parton distribution functions |], §, || || . The latter possess hybrid properties: in different regions 
of phase space they share features with ordinary parton densities of deep inelastic scattering 
and with exclusive distribution amplitudes. They smoothly interpolate between these two limits. 
At the present stage all the required perturbative inputs are available, i.e. one-loop coefficient 
functions |5j, [| [7], [8|, |§ and two-loop anomalous dimensions of the moments of the non-forward 
functions [[10], |TTJ, which make possible its analysis in next-to-leading order (NLO) of perturbation 
theory in the flavour singlet channel. In the above list of corrections the latter were derived within 
a formalism which allows to determine the eigenf unctions of the two-loop generalized Efremov- 
Radyushkin-Brodsky-Lepage (ER-BL) evolution equations in closed analytical form. This is a 
crucial result which is not accessible in the direct calculation of NLO kernels. 

The leading order solution of the evolution equations, which are governed by the kernels 
evaluated in Refs. 0, |^, O [I3L Rl [IS] , in terms of the conformal partial wave expansion was 
given in Refs. f||, |^, |l(| (see also [|17| and for a direct numerical integration of evolution equation 



Ref. QT2j] ) . In two- loop approximation the corrections to the eigenf unctions were derived in || |TT 
In the present investigation we will use these results for an explicit study of the magnitude of 
these effects on the evolution of the flavour non-singlet non-forward parton distribution function 
versus the LO evolution considered by us previously |R) 



2. Solution of the two-loop evolution equation. To start with let us first describe the formalism and 
spell out our conventions. Here we adopt the definition of the non-forward parton distributions 
introduced by Radyushkin |2|, |3| which fulfill the evolution equation 

/i 2 ^0(z, 0=J dx'K (x, x', C \a s (Q 2 )) 0(at, (), (1) 

where the kernel K (x, x 1 , ( \a s (Q 2 )) can be calculated order by order in perturbation theory. 
Since the leading order evolution equation can be diagonalized with the help of the conformal 
operators it is convenient to employ this partial conformal wave expansion also beyond leading 
order although the Gegenbauer polynomials are not the eigenfunctions of the two-loop generalized 
ER-BL equation. Namely, the solution of the Eq. (|l|) can be written in the form 

oo 

0(a;,C,Q 2 ) = E^(^C|« s (g 2 )) OjiCQ 2 ), (2) 
where the multiplicative renormalizable moments evolve as follows 

OACQ 2 ) = exp {-- £f (a s (r))\ O^CQl), (3) 
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with the forward anomalous dimensions we need at 0(a 2 s ) accuracy 7° (a s ) = ^ 7 J -°' ) + jj 1 ^ + 
. . ., where = -C F (3 + (i+1) 2 (i+2) - + 2) + 4^(1)) and 7^ can be found in Ref. @. We 
have chosen the initial condition so that there are no radiative corrections at the low normalization 
point Ql so that the Oj((,Ql) are given by ordinary Gegenbauer moments of the non-forward 
distribution which are related to the matrix elements of the tree level conformal operators by 



dx cf (2^-1] 0(x, c, Ql) = ^(ti\ ,-(}<). yi-cy 1 { D + /d + ) v> 



X 



1 



3/2 



c 



\h). (4) 



The problem is reduced to finding the correction to the eigenfunction. It was solved in our previous 



studies U M, 11, 19 



Z7T k=j+2 

where the LO partial conformal waves are defined via the Gegenbauer polynomials 



Ixl x 



cT 2--1 



with Nj 



U + W + 2) 
4(2j+3) 



(5) 
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The function <P jk is §, [H| y, [19 



$jk (a s (Q 2 )) = Sjk (« s (g 2 )) {d jk ( 7 £ 0) - ft) - <7; fe } 



(7) 



with 



djk 
9jk 
with A 



-V(-iy- fe ] 

2 [ +l J J 0'-^)(j' + ^ + 3)' 
2C F 1 2A jfc + (A,-* - + 2) + V(l)) 



C/-A;)0" + fc + 3) 



if, 



'j + k + A" 



(k + l)(k + 2) J ' 
+ 2V0'-A;)-V0' + 2)-V(l)- 



(9) 



The factor Sjk appears as result of the evolution of the coupling constant and reads 



Sjk («s(Q 2 



(0) (0) 

Tj -Ik 

(0) (0) . n 



1 - 



MQ 2 ). 



2)Xl+(7f-7f)^o\ 



(10) 



with A) = %T F N f - fC A . 

Obviously, taken as they stand the above Eqs. (@,H,|D are valid only for the non-forward 
distributions with support < x < ( since the Gegenbauer polynomials Cj(2x — 1) form a 
complete set only on the interval x e [0, 1]. But as has been noted in Ref. || that means that we 
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can understand the above expansion only in a restricted sense. Namely, we have to represented it 
in the form HI 



x f x 



U ~ 1/2 C» (2- 2 1 - 2 - r (l) r 0' + 2 ^) f 1 dtm ^-i/2 5 ( j)(ct _ 



and treat the RHS as a mathematical distribution^. In order to circumvent this disadvantage 
we expand the non-forward distribution in a series of polynomials Vj(x) orthogonal in the region 
< x < 1. Then the expansion coefficients will be given as convolution of Vj(x) with the RHS of 
Eq. (0). 

Note that we can exploit any set of orthogonal polynomials Vj(x) available. Namely, the 
general expansion looks like 

00 

<D(x,(,Q 2 ) =£^(x)A^(C,Q 2 ), (12) 

3=0 

where the conjugated polynomials Vj(x) are defined such that 

-1 

dx Vj(x)V k (x) = S jk . (13) 
And the moments A4j(£, Q 2 ) are given by a finite sum 

Mj((,Q 2 ) = £Ef k (OO k ((,Q 2 ), (14) 

k=0 

with Oj-moments expressed in terms of the original ones as follows 

Ojic, q 2 ) = a,(c, q 2 ) + E % ( a s(Q 2 )) o*(c, q 2 )- (15) 

Zir k=0 

Taking into account the above observation, the expansion coefficients are given by the integral 

Ef k (0 = C dx ^C 3 k /2 (2x - l)V 3 ((x). (16) 
Jo l\ k 

Let us repeat that we can use any appropriate polynomials for this purpose. The criterion for 
choosing a specific one is thus the fastest convergence of the series. Below we give the result for 
the expansion coefficients of the Jacobi polynomials. Namely, 

z a x^ 
nj(a,py j 



Vj(x) = P^{2x - 1), Vj{x) = X . X Q . P^\2x - 1) 



with nAa,p) = -. , 7=-. r 17 



2 For this reason we omit the spectral constraint which appears on the LHS as a result of integration. 



The expansion coefficients can easily be obtained by the methods developed in Ref. [IT]] and read 

T(A; + 2) T(j + P+1) T(j + k + a + p + l) 




3) T(k + (3 + 1) T(j - k)T(j + a + (3 + 1 
■j + k,j + k + a + p + l,k + 2 
2k + A,k + P + l 



C ■ (" 



The results for all other classic orthogonal polynomials immediately follow from this expression. 
For special values of the parameters the Jacobi polynomials coincide [^0| either with Gegenbaueij], 
Pj X 2 ' X 2 \x) = r^^rxA+i) ; or LegendreQ, Pj 0,0 \x) = Pj(x), or Chebyshev polynomials 

of the first, pj~*'~*\x) = ^^-T^x) and second kind, P^\x) = ^^L Uj{x). 

The solution of the renormalization group equation (|3]) in two-loop approximation can be 
written in the form 




I Po + Pi 4n \ ' 



JO) 



Ij I A) ( la JO) Jl)\ n (rt 2\ n (ra 



•Itt 



/3 l7 f it \ oc s {Q 2 )-a s {Q 



2P 2 Po / 2tt 



(19) 



where the expansion in the second line is done in order to treat two-loop corrections to the 
evolution on the same footing as one-loop corrections to the Wilson coefficients, which when both 
are summed in the amplitude allows to minimize the renormalization scheme dependence (see, 



for instance, the first paper in Ref. [18] and [21]). The coupling constant in NLO of perturbation 
theory can be approximated by 

where (5\ is the second coefficient in the expansion of the QCD /^-function | = ^ Po+(t*) + - ■ ■ 
and it reads ft = fC A N f + 2C F N f - f C\. 

3. NLO evolution of the model distributions. In this section we will use the results given above 
for explicit numerical studies of the evolution of the non-forward parton distributions. To this 
end we need an initial condition for the evolution equation. The most adequate model for the low 



scale input functions was given by Radyushkin in Ref. ||22j| . Namely, 0(x,Q is defined in terms 



of the double distribution function F(y, z) H via the following relation 

0(x, C, Ql) = I' dy f 1 dz F(y, z, Q 2 ) 9(1 - y - z)S(x -y-(z). (21) 
Jo Jo 



3 We have use them in our preliminary study of LO evolution effects 16 . 

4 This possibility has been discussed in Ref. |]l7| but the authors did not manage to find an explicit analytical 
expression for the expansion coefficients. 
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Here we have omitted the dependence on the t-channel momentum transferred squared. Although 
for a massive target the formal limit t — ^ is not accessible for nonvanishing £ due to the 
kinematical restriction, m|£ 2 /£ — — t [23, ||, this condition does not affect the results of the 
evolution. 

For our study we accept the following model for the double distribution function corresponding 
to the non-singlet function 0(x, £, Ql) |p2| 

F(y,z,Ql)=q(y,Q 2 )7c(y,z), (22) 

with the plausible profile functions ir(y, z) 

n(y,z)=6^(y-z). (23) 

In Eq. (^21) the function q(y, Qq) is an ordinary forward parton density measured in deep inelastic 
scattering taken at a low normalization point. We will also consider the asymptotic distribution 
functions which although phenomenologically probably irrelevant serves as a good probe for the 
net evolution effects in NLO approximation since its scale dependence is governed by the non- 
diagonal elements of the anomalous dimension matrix of the conformal operators which appear 
only beyond leading order. It is given by the first term in the expansion (^) |2"2"| , IE| 



O u (x, C, Ql) = 6^1 - 9(( - x)O ((, Ql), (24) 

(see Fig. [I] (a)) where Oo(CQo) ( see ^q- dD) * s a m& trix element of the conserved local (axial- 
hector current, 0o(C>Qo) = (^ / |^(0)-^' ? /'(0)|^) with r = 7+, 7+75. Therefore, it depends neither 
on the skewedness parameter nor on the renormalization point Qq. 

Since in present paper we are studying only the non-singlet evolution we consider the combi- 
nation 

q NS (x) =u{x) -d{x), (25) 

and take as initial condition the following CTEQ4M parametrizations for the u and <i-quark 
distributions P|| 



u(x) = 1.344x- a499 (l-x) 3 - 689 (l + 6.042a; a873 ), (26) 
d(x) = 0.640x- a499 (l-x) 4 - 247 (l + 2.690x°- 333 ), (27) 

at an input scale Qo = 1.6 GeV. Then the above model reads form these and is given in Fig. ^| 
(a). Note that the factor Oq((, Ql) in the asymptotic distribution is defined by the first moments 
of q NS (x), Co (C, Ql) = Jo dx q NS (x). Let us mention that the models we use satisfy the positivity 
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Figure 1: Evolution of the asymptotic distribution function. The x — (^-dependence given by Eq. 
fl24|) is shown in (a). The input distribution (dashed-dotted line) was evolved in NLO approxima- 
tion (solid line) up to Q 2 = 100 GeV 2 for the skewdness parameters £ = 0.1 (b), ( = 0.5 (c) and 

C = i.o (d). 



constraints^ derived in Ref. |2^, |22 



O(x,0 < sjq{x)q ((x - C)/C) Tl 



f28) 



Moreover, it is useful to mention that Eq. (|2T|) with profile defined by (|23|) saturates the constraint 
inequality ( p8|) in the region of its validity x > (. 

Now we are in a position to present the numerical results for the evolution of the models 
introduced above. The appropriate values for the parameters which we did not mention so far 
are Nf = 4, Af^g = 220 MeV. For the NLO anomalous dimensions we have used the simplified 



expression derived by Yndurain et al. [[18] which works with an accuracy better then 0.2%, namely 
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(i) 



A 2-^1 



(3 + I) 



(29) 



5 Note an extra factor of l/\/C i n Eq. ( ^ ) found in Ref. (22) which was missed in p5f . 
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Figure 2: 3D shape of the Radyushkin's model (a). Evolution of input distribution (dashed- 
dotted line) in leading (dashed lines) and next-to- leading order (solid lines) with ( = 0.1 (b), 
C = 0.5 (c), C = 1.0 (d) and Q 2 = 10 GeV 2 (curve a), Q 2 = 100 GeV 2 (curve b). 
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In the following we will have a closer look on the distribution functions evolved up to the 
reference scales Q 2 = 10, 100 GeV 2 for the skewedness parameters set equal to ( = 0.1, 0.5, 1.0 and 
compare them with LO results. We exploit for these purposes an expansion in terms of Legendre 
polynomials. We perform the evolution^ by evaluating 70 moments in the series (|T2"| ) in the case 



3 The calculations were done with a code written for MAPLE. 



of the Radyushkin's model distribution. In the particular case of the asymptotic distribution the 
calculation was made in a different way. Due to the support properties of the latter we used 
the expansion of the non-forward distribution in Gegenbauer polynomials C^ 2 (2x/( — 1), i.e. we 
directly employed Eqs. (fj,^,^]). Thus there is no need for a double expansion which would restrict 
the number of terms one can treat in the expansion. We have used up to 100 polynomials. In 
both cases we made a fit to get smooth curves instead of rapidly oscillating ones. 

The results are shown in Fig. |]||. As can be seen there is only a very small difference between 
NLO and LO evolved distribution. As expected this deviation grows with increasing Q 2 , but it 
remains small even for large Q 2 . Note that for asymptotically large Q 2 the excitation of higher 
harmonics will die out and both distributions take the form O as (x, () from Eq. (f2T|). 
4- Conclusion. To conclude, in this note we presented our results on the numerical evolution of the 
non-forward distribution function in two-loop approximation in the flavour non-singlet channel. 
We studied two models: the asymptotic function and Radyushkin's model f2"2]| , we have found 
that the net effects of NLO corrections to the evolution kernels are extremely small and do not 
exceed the level of a few percent for moderately large Q 2 . These results together with the fact 
that the evolution of the (9 as (:c, () (|24"D is governed by off-diagonal elements of the kernel in the 
basis of Gegenbauer polynomials suggest that the latter are small as compared to diagonal entries 
and, therefore, the kernel is quasi- diagonal in this basis. 
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